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Abstract 

We study the four dimensional Gaussian spin glass in presence of a magnetic 
field. Using off-equilibrium numerical simulations we have found that the probability 
distribution of the overlaps is built in the same way as that of the Mean Field ap- 
proximation with replica symmetry breaking. Finally we have studied the violation 
of the fluctuation-dissipation theorem in presence of magnetic field. 



1 Introduction 



One of the main open questions in the study of finite dimensional spin glasses is the 
existence of a phase transition in presence of an external field. While in absence of magnetic 
field the droplet model Q and the Mean Field (MF) picture || both predict the existence 
of a phase transition, in presence of magnetic field the situation is fully different: Mean 
Field predicts a phase transition whereas the droplet model shows that the magnetic field 
destroys the frozen phase. 

Moreover there are some analytical arguments which imply that the phase transition in 
presence of magnetic field is of a rather peculiar type. For instance Bray and Roberts M, 
working with a reduced theory (obtained by projecting the original theory in the replicon 
subspace; in presence of magnetic field that is the only critical mode), have shown that there 
is no weak-coupling fixed point in magnetic field in their renormalization group equations 
also near 6 dimensions. An absence of a weak-coupling fixed point is often taken as an 
indication of a first order transition. Here the situation is quite less clear. 

The very existence of a transition is still controversial and not too much work has been 
devoted to its study. Numerical simulations have been done in the past [Q. They were 
compatible with the possibility of a transition, but the situation was not so clearly cut and 
no convincing conclusions could be reached. Only recently an off-equilibrium numerical 
simulations in 4 dimensions have strongly suggested the existence of a phase transition in 
presence of magnetic field ||. 

Two main advantages in using dynamical methods are that we can simulate very large 
systems (up to 40 4 in this work) loosing practically all the finite size effects and that they 
are quicker with respect to an equilibrium simulation, because we don't need to thermalize. 
These methods have been already largely used in the numerical studies of spin glasses (see 
for instance reference ||). 

In the broken replica symmetry solution of the SK model || 0] there is a phase transition 
also when the system is plunged into a magnetic field and the line T c (h) which separates 
the paramagnetic from the spin glass phase is called Almeida-Thouless (AT) line ||. The 
order parameter of the Mean Field theory is the probability distribution of the overlaps, 

In absence of external field P(q) is a delta function centered on q = for T > T c , while 
for T < T c it becomes a highly non trivial function with two delta functions on the values 
q = ±^ea (<?ea is the maximum allowed value for the overlap) and a non-zero part between 
them. 

When a magnetic field is switched on the function P(q) becomes zero for every negative 
overlap and the minimum allowed value for q is shifted upwards from — ^ea to q m m > 0, 
while the maximum value (^ea) almost doesn't change. This means, in terms of the 
distribution function of the overlaps, that for T > T c (h) the function P(q) is a delta 
function centered on a strictly positive value and that for T < T c (h) the function P(q) is 
the sum of a delta function on the maximum value g max = qea plus a non-zero part down 
to g m i n > and a smaller weighted delta function on g min . 

In this paper we present evidences for a Mean Field like phase transition at finite 
temperature. We show that the order parameter P(q) has a non-zero support: we have 
characterized numerically the mean, the minimum and the maximum value allowed, de- 
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noted q, q m i n and g max respectively, and we have found that g m j n <q< q max . 

The plan of the paper is the following. In the next section we fix the notation and we 
describe the quantities we have measured. In sections three we show the numerical results. 
Finally we present the conclusions. 

2 The model and the observables 

We have simulated the Gaussian Ising spin glass in four dimensions on a hypercubic lattice 
of volume N = L A with periodic boundary conditions. The Hamiltonian of the system is 
given by 

Ti — - ^ OiJijOj - ^Y^oihi . (1) 

<ij> i 

By < ij > we denote the sum over nearest neighbor pairs. The are Gaussian variables 
with zero mean and unitary variance. The external field is also Gaussian with zero mean 
and variance h\. We have studied systems with ho = 0.05, 0.1, 0.2, 0.3, 0.5. 

We can justify the choice of a Gaussian magnetic field as follows. The starting point is 
the Ising spin glass Hamiltonian with an uniform magnetic field h 

H = - Y (TiJijO-j - h ^2<Ji . (2) 

<ij> i 

We can perform the following transformation (a "local gauge transformation") on the 
couplings: J i3 - — > J-j = UiJ^n^ where rii = +1 or —1. This transformation leaves the 
Hamiltonian 7Y unchanged because the probability distribution of the couplings is Gaus- 
sian (Jf- = (J'ij) 2 ). Now we recast the spins to = nj<7j, and finally our Hamiltonian 
reads 

7-Lo = ~ s i J ij s j ~ Yl( h o n 'i) s i » ( 3 ) 

<ij> i 

where n-rij = 1. We remark the full arbitrariness of the choice of n-'s. In particular we can 
choose them from a bimodal distribution: i.e. n[ = 1 with probability 1/2 and —1 with 
the same probability. And so, if we define h[ = hon'^ the Hamiltonian can be written as 

Ho = - J2 s i J ij s 3 - J2 h 'i s i ■ ( 4 ) 

<ij> i 

We have therefore shown that a spin glass with an uniform magnetic field (ho) is equivalent 
to a spin glass in which the magnetic field is random with zero mean and variance h^. The 
probability distribution of such a magnetic field is bimodal, not Gaussian. Nevertheless 
there are not reasons to suppose a different physical behavior of these two cases (bimodal 
and Gaussian). We have chosen a Gaussian distribution, and not a bimodal one, because 
in the Gaussian case there are exact relations among some quantities. 

We are interested in measuring the mean value of the overlap between two replicas, 
which is defined as q = J qP(q)dq, without doing expensive equilibrium simulations. If 
we take two replicas in random configurations (as we do at the beggining of a simulation) 
their overlap is zero (q(t = 0) = 0). Letting them evolve, the overlap will never increase 
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beyond q m i n , denned as the minimum overlap allowed at the equilibrium []. This fact has 
been largely verified also in off-equilibrium simulations without magnetic field: during the 
simulation the overlap fluctuates around zero or slightly grows. This observation gives us 
a practical tool to calculate g m i n via an infinite-time extrapolation, but also asserts that we 
cannot get information on the whole P(q) simply by looking at the off-equilibrium overlap. 

To measure q we have exploited a relation valid at equilibrium when the applied field 
is Gaussian, which reads 

((Tihi) 1-q . 

~w = — ' (5) 

where with the overline we mean an average over the quenched random interactions and 
external fields. This relation can be easily obtained via an integration by parts (exploiting 
that hi is a Gaussian random variable) and it's exact also in a finite volume. In order 
to compute q we can measure (<Ti(t)hi), which is a quantity that rapidly converges to its 
infinite-time value. The fact that g min differs from q is a clear signal of replica symmetry 
breaking. 

The second part of our study is focused on the fluctuation-dissipation theorem (FDT) 



and its generalization in the out of equilibrium regime [10, O, [12|, called off-equilibrium 



fluctuation-dissipation relation (OFDR). In reference flT2[ a detailed study of such a relation 



in finite dimensional spin glasses without magnetic field can be found. Here we extend those 
studies in presence of a magnetic field, obtaining similar results, which confirm the Mean 
Field behavior of the phase transition. 

To study the OFDR we have measured the spin-spin autocorrelation function, C(t, t w ), 

and the integrated response of the system, (aihi)/el, where the perturbation to the Hamil- 
tonian 7i (eq.([l|)), h, is a random Gaussian magnetic field with zero mean and variance €q. 
In the next paragraphs we will obtain a formula that links, even in the early times of the 
dynamics, the response and the auto-correlation function. 

Given a quantity A(t) that depends on the local variables of our original Hamiltonian 
(7i). We can define the associate autocorrelation function 

C(t,t') = {A(t)A(t')) , (6) 

and the response function 

(7) 



6e(t>) 

where we have assumed that the original Hamiltonian has been perturbed by a term 



W = H + J e(t)A(t) dt . 



The brackets ((■■•)) ineq.(§) and eq.(|7|) imply here a double average, one over the dynamical 
process and a second over the disorder. 

lr This can be better understood if we think that P(q) is the equilibrium distribution of a dynamical 
variable which feels a potential V(q) = — log P(q). The peak at q = q m i n in P(q) corresponds to a deep 
well in V(q) and if q starts with a value less than g m j n it will go downhill towards g m i n - See reference 0] 
for more details. 
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In the dynamical framework assuming time translational invariance it is possible to 
derive the fluctuation-dissipation theorem, that reads 



R{t,t')=m-t')^P- ■ (9) 

In spin models a common choice for A(t) is A(t) = c?i{t) or A(t) = iV _1//2 Y,i Oiif)- in 
this case, because the system feels a magnetic field, to have a simpler response we should 
perturb it with a random field, h, and measure the staggered magnetization. In order to 
derive a fluctuation theorem where the response is related to the one site correlation, we 
must chose a perturbation such that the off-site elements of the response are zero 0. So here 
we put A(t) = N~ l l 2 Y.i ni&i(t), where rij = hi/e . Thanks to the fact that (niUj) = 5 it j, 
we have that with this choice 

c(t,t') = W(^(t)^(t')) (10) 

i 

and 
where 

rn s [h](t) = ^Y,(nMt)) (12) 

i 

is the staggered magnetization, which is a functional of the magnetic field, h(t), and a 
function of the time. 

The fluctuation-dissipation theorem holds in the equilibrium regime, but in the early 
times of the dynamics we expect a breakdown of its validity. Mean Field studies fll"3| 
suggest the following modification of the FDT: 

R(t,t') = (3X{t,t')9{t-t')^^ . (13) 

It has also been suggested in fl~3| , [14|, |15| that the function X(t,t') is only a function of 
the autocorrelation: X(t,f) = X(C(t,t')). We can then write the following generalization 
of FDT, which should hold in early times of the dynamics, the off-equilibrium fluctuation- 
dissipation relation (OFDR), that reads 

R(t,t') =px(C(t,t'))8(t-t')^^ . (14) 



We can use the previous formula, eq. ([14]) , to relate the observable quantities defined in 
eq.flTUD and eq.(|T^). Using the functional Taylor expansion we can write 



/CO 
dt 
-co 



, 5rn s [h)(t) 



5h(t>) 



h(t') + Oih 2 ). (15) 

h(t)=0 



2 We notice also that in absence of a magnetic fields, a constant field perturbation is gauge equivalent 
to a random magnetic field. 
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Exploiting the definition of eq. flTT|) and using the fact that m[0](i) = for every perturbing 
field orthogonal to the pre-existing one, i.e. such that Y^,i( n ihi) = (which is true if hi is 
another (uncorrelated) random field, as happens in our case), we obtain 

m s [h](t) = f dt' R(t,t')h(t') + Oih 2 ) . (16) 

This is just the linear- response theorem neglecting higher orders in h. 

By applying the OFDR we obtain the dependence of the staggered magnetization with 
time in a generic time- dependent magnetic field (with a small strength), /i(t)|] 

m s [h](t) * P f jt' X[C(t,t')]^^-h(t') . (17) 

Now we can perform the following experiment. We let the system evolve with the unper- 
turbed Hamiltonian of eq.(P from t = to t = t w , and then we turn on the perturbing 
magnetic field h, which is Gaussian distributed with zero mean and time-independent 
variance, eg. Finally, with this choice of the magnetic field, we can write?] 

m a [h](t) - *o0jf ^ X[C(t,t')]^p- , (18) 
and by performing the change of variables u = C(t,t'), equation (|T8|) reads 

m s [h){t) ~ e p [ duX[u), (19) 

JC{t,t w ) 

where we have used the fact that C(t, t) = 1 (always true for Ising spins). In the equilibrium 
regime (FDT holds, X = 1) we must obtain 

m s [h)(t) ~e p(l-C(t,t w )) , (20) 

i.e. m s [h](t) T/eo is a linear function of C(t,t w ) with slope -1. We remark that we can use 
this formula to obtain g max as the point where the curve m s [h](t) versus C(t, t w ) leaves the 
line with slope — /3e (as we will explain). 

In the limit t, t w — * oo with C(t,t w ) = q, one has that X{C) — > x(q), where x(q) is 
given by 

x(q) = I' dq' P(q') , (21) 

where P(q) is the equilibrium probability distribution of the overlap. Obviously x(q) is 
equal to 1 for all q > q max , and we recover FDT for C(t,t w ) > q max . This link between 
the dynamical function X(C) and the static one x(q) has been already verified for finite 
dimensional spin glasses ||12|| . 

For future convenience, we define 

S(C) = I' dq x(q) = f 1 dq f dq' P(q') . (22) 

JC JC Jq min 

3 The symbol ~ means that the equation is valid in the region where linear-response holds. 
4 We ignore in our notation the fact that m s [h] (t) depends also on t w . 
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or equivalently 



In the limit where X 



P(q) 



d 2 S(C) 



d 2 C 



C=q 



x we can write eq.(|T9"D as 
m s [h}(t)T 



S(C(t,t w )) . 



(23) 



(24) 



Looking at the relation between the correlation function and the integrated response 
function for large t w we can thus obtain q maiX , the maximum overlap with non-zero P(q), 
as the point where the function S(C) becomes different from the function 1 — C, and g m i n 
as the smallest value of C. 

At this point we have numerical methods to compute three important different values 
of q: g ma x, <?min and q. 



3 Numerical results 

3.1 q and q m[n 

We are interested in the behavior of the system in the out of equilibrium regime, so we do 
not need to thermalize the sample and we can simulate very large samples of millions of 
spins (24 4 , 32 4 and 40 4 ). We expect our data not to be affected by large finite size bias: 
we find that, in the range of temperature considered, the data for different lattice sizes 
(L = 24, 32, 40) coincide within the error-bars, with the largest systems (L = 32, 40) giving 
practically the same values. 

All the numerical simulations have been performed on the parallel super-computer 
APE100 [T§. 



In the first part of our study we have done simulations using the Hamiltonian of eq.(JJ) 
without perturbing the system (h = 0). 

According to our dynamical approach, we are interested in doing measurements in the 
off-equilibrium regime of large times. So we perform the simulations following an annealing 
schedule, with slower and slower cooling rates. This is equivalent to run a simulation at a 



fixed temperature, considering larger and larger waiting times |T2j , PH . The advantage is 
twice: within a single run we are able to collect data at different temperatures and there are 
smaller finite-time effects because when the temperature is lowered by an amount AT the 
time needed to "forget" the previous temperature is smaller than if it had been quenched 
from T = oo. In the present case the range of temperatures [] is from T = 3.0 down to 
T = 0.5 with a step of AT = 0.25. The cooling rates have been chosen in such a way 
that in the s-th annealing run the number of Monte Carlo Steps (MCS) is proportional to 
2 s , with s ranging from to 12. The annealing procedure is the following: indexing the 
temperatures from the highest (zt=3.o = 1) to the lowest («t=o.5 = H) the number of MCS 
at each temperature is 2 s • i\. This means that in the slowest run the system will stay at 
the lowest temperature for about half a million MCS. 

5 We remind that in the four dimensional Gaussian Ising spin glass without external field the transition 
temperature is T c = 1.8 [p7| . 
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Figure 1: Data of q(t) (bottom) and 1 — (ai(t)hi) T jh\ (top) versus Monte Carlo time. 
The curves are the best fits (see text). L = 40, h = 0.3 and T = 0.5. 



To analyze the data, we fix a temperature, we recollect the data at that temperature 
from different annealing runs and we extrapolate the result in the limit of infinitely slow 
cooling. This will give us information on the large times off-equilibrium regime. 

In particular we look at the overlap between two replicas, q(t), which start very far 
from each other: q(t = 0) = 0. In the limit t — > oo this overlap will tend to the minimum 
overlap allowed at the equilibrium, q m i n . Also we have measured ((Ti(t)hi) from which we 
obtain the value of q using eq.(Q). In figure [I] we show the data (L = 40, h = 0.3 and 
T = 0.5) whose infinite-time extrapolation gives the mean and the minimum value for the 
overlap (top and bottom data respectively). 

We fit the data plotted in figure [1] with the following formula: q(t) = At~ B + C. If 
the value of the exponent B were too small the uncertainty on the value of C should be 
very large and the significance of the fit very poor. We have found that it decreases with 
decreasing temperature or magnetic field and we have checked in all our fits that it were not 
too small. In figure [I] we present the worst case we can fit satisfactory (lowest temperature 
and ho = 0.3), obtaining as the best parameters: A = —0.069(2), B = 0.20(3) and 
C = q = 0.792(5) for the top data and A = -1.4(2), B = 0.20(3) and C = q mill = 0.64(2) 
for the bottom data. Both fits have very good \ 2 values. Note that the best B exponent 
is the same in both fits. The presence of a power law approach to equilibrium can be 
associated to the existence of a well defined off-equilibrium correlation length in finite 



dimensional spin glasses which grows with a power law of the time [0, [19 



The dynamical approach in presence of a magnetic field has the advantage that, using 
a fitting function of the type At~ B + C for any observable in the off-equilibrium regime, the 
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Figure 2: g (top) and g min (bottom) versus T for L = 40 and ho = 0.3. 

best value for the exponent B is always greater than in the case of exactly zero external 
field and moreover it has a finite limit when T — > 0, while for h = very often B oc T . 
The reasons for this difference in behaviour are unclear to us Q 

Now we clearly see the usefulness of the trick of using a Gaussian magnetic field: we 
can calculate q from the extrapolation of (<Ji(t)hi), which is a practically constant quantity 
with respect to q{t) (look at figure [I]). Moreover we don't need to do any limit of small h 
and so we don't have the problem of very small B exponents. 

The same kind of fitting analysis has been done for all temperatures, obtaining values 
for g m i n (T) and q(T), plotted in figure ||. In the high temperature phase the extrapolated 
values for g min and q coincide (as they should because the equilibrium overlap distribution 
function is a delta function) confirming the correctness of the dynamical method. 

Figure |2] gives clear evidence of a wide region (T < T c (h = 0.3) ~ 1.2) where the order 
parameter P(q) is not a single delta function. We have calculated the functions q m m{T) 
and q(T) also for different values of the magnetic field. We see a clear bifurcation also 
for ho = 0.5 at a temperature T c (h = 0.5) ~ 1.0, while for lower magnetic fields the 
results are less clean because the errors are larger (due to the just described problem in 
the extrapolation procedure). 

The last result we present in this section is the shape of the equilibrium magnetization 
in a magnetic field as a function of the temperature. More precisely what we show in fig |^ 

6 We would like to recall that numerical simulations of relative small size systems for the SK model 
indicate that the time for passing at equilibrium from one value of q to an other value of q with the same 
sign increases as exp^iV 1 / 4 ) while that for changing the sign of q increases as exp(AN 1 ^ 2 ). In magnetic 
field only the first processes are present and in the SK model they correspond to much faster movements. 
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Figure 3: (aihi)/h% versus T for L = 40 and h = 0.3. The curve is only a guide to the 
eye. 



is the staggered magnetization {aih^) jh\^ but we have verified that it behaves qualitatively 
and also almost quantitatively like the magnetization in a constant field. 

We observe a behavior similar to the one measured in experiments on real samples: the 
magnetization has a little peak, whose height is greater than the zero-temperature value 
by a few percents. We want to stress that the peak is at a temperature higher than the 
critical one, so the magnetization peak is not exactly on the AT-line. 

3.2 Fluctuation-Dissipation in presence of magnetic field 

In the second part of our study we have simulated an Ising spin glass in a Gaussian mag- 
netic field with variance (ho = 0.2,0.3,0.5) at a fixed temperature (T = 0.75, 1.0, 1.5). 
We have measured, for various waiting times (t w = 2 8 , 2 11 , 2 14 , 2 17 ), the autocorrelation 
function, defined in eq. (|10D , and the staggered magnetization, defined in eq. (|12"D , with dif- 
ferent amplitudes of the perturbing field (eo = 0.02,0.03,0.06). We are interested in the 
relation between these two quantities, which tends to the function X(C) in the large t w 
limit (see eq . (|24"D ) . The data do not depend on the choice of eo or L, so the results we will 
show are in the linear-response regime and they are not affected by large finite size biases. 
In these runs the starting configuration of the spins is always random. 

We can keep in mind the possible link between the static and dynamics (that we do 
not study in the present paper) for a spin glass in presence of a magnetic field and we 
can compute the qualitative form of the function X(C) using as input the S{C) function 
obtained in the Mean Field approximation and in the droplet model. 
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Figure 4: m s [h](t) T/e versus C(t,t w ), with T = 0.75, h = 0.3 and e = 0.06. The 
straight line is 1 — C . 

From the function S(C) we can get information on the overlap distribution function 
P{q), through eq.([23|). Let us remind which is the prediction for the S(C) assuming the 
validity of one of the competing theories described in the introduction. The droplet model 
predicts P(q) = S(q — q) and consequently 

w={±: q c tail: ^ 

i.e. there is no dependence of the staggered magnetization on C(t, t w ) in the off-equilibrium 
regime (C < q), like an ordered ferromagnet PD[ . On the other hand the MF like prediction 
for the overlap distribution P(q) = (1 — XM)8(q — g max ) + x m 5(q — q m i n ) + p(q) (where 
the support of p(q) belongs to the interval [g m i n , <?max]), implies that 

f S(0) for C < q min , 
S(C) = I 5(C) for q min <C< g max , (26) 
I 1-C for C > g max , 

where s(C) is a quite smooth and monotonically decreasing function such that 

C=q 

In figure | we plot m s [h](t) T/e versus C(t,t w ), with T = 0.75, h = 0.3 and e = 0.06. 
The data are the average over 6 samples of a 32 4 system. We can see from figure § that, 
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Figure 5: m s [h](t) T/e versus C(t,t w ), with T = 1.5, h = 0.3 and e = 0.02. The straight 
line is 1 — C. See text for the short vertical and horizontal lines. 

when the data are not on the straight line (FDT or quasi-equilibrium regime), they do not 
lie on an horizontal line, i.e. they depend on the value of the autocorrelation also in the 
off-equilibrium regime. So we can conclude that the droplet theory is not able to describe 
the data in the frozen phase. 

Interpreting the data using the MF picture, we deduce from the small curvature of the 
function S(C) in the region C < g max [] that p(q) is small. Anyway the existence of two 
different delta functions in P(q) is clear. 

To verify the correctness of the method we present in figure |5] the same kind of plot for 
a system which is in the paramagnetic phase (T = 1.5, h = 0.3 and eo = 0.02). The data 
stay, as they should, on the equilibrium line S(C) = 1 — C. The few points which leave 
the straight line are those with the lowest t w and large t. This effect is due to the fact 
that the autocorrelation C(t,t w ) tends to q m i n only if the system at time t w has reached 
equilibrium. On the contrary if the system is in a random configuration at time t w the 
autocorrelation will tend to zero also in a magnetic field. In our simulation the system is 
in an intermediate situation and so the lowest value for the autocorrelation with t w = 2 11 
is something smaller than g m i n , while for t w = 2 17 we think that the system is nearly 
equilibrated and the autocorrelation does not decrease beyond q min . In figure [| we have 
reported the value for the equilibrium overlap, q = g min = g max ~ 0.356, calculated with the 
annealing runs (vertical line) and the corresponding staggered magnetization (horizontal 
line). Though it is not very clear from figure [5], we have verified that for every t w the 
staggered magnetization saturates to the value marked with an horizontal line and that 

The data can be fitted also with a straight line, whose non-zero inclination is almost i^-independent. 
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for the greater waiting time the autocorrelation tends to the vertical line. 



4 Conclusions 

We have obtained, using off-equilibrium simulations, g m j n , g max and q for the four dimen- 
sional Gaussian spin glass in presence of a magnetic field finding that in the low temperature 
region g min <q< g max according with the predictions of Mean Field Theory. 

This result point clearly toward to a phase transition between a spin glass phase with 
spontaneously broken replica symmetry (g m i n < q < g max ) and a phase where the replica 
symmetry is stable (g min = q = g max ). 

Moreover we have extended the numerical studies of the validity of the fluctuation- 
dissipation theory obtaining that the function that determines the violation depends only 
on the correlation, as Mean Field theory predicts. We plan in the future to try to link this 
function (X) with that obtained from the static of the system (x) as it has been done in 
absence of magnetic field although some indirect evidences of this link has been shown in 
the present paper. 
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